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Abstract 



Using the renormalization group method, we improved the first order solution 
of the long- wavelength expansion of the Einstein equation. By assuming that 
the renormalization group transformation has the property of Lie group, we 
can regularize the secular divergence caused by the spatial gradient terms and 
absorb it to the background seed metric. The solution of the renormalization 
group equation shows that the renormalized metric describes the behavior of 
gravitational collapse in the expanding universe qualitatively well. 
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I. INTRODUCTION 



The spatial gradient expansion of the Einstein equation is non-hnear approximation 
method which describes long-wavelength inhomogeneity in the universe. This approxima- 
tion scheme expands the Einstein equation by the order of the spatial gradient. As the 
background solution, we solve the Einstein equation by neglecting all spatial gradient terms. 
The resultant solution has the same form as that for the spatially flat Friedman-Robertson- 
Walker(FRW) universe, but the three metric can have spatial dependence. It is possible to 
include the effect of spatial gradient terms by calculating the next order. This method can 
describes the long-wavelength non-linear perturbation without imposing any symmetry for 
a space, and is suitable for analyzing the global structure of an inhomogeneous universe. 
Furthermore, it is easy to include perfect fluid with pressure and scalar fields. 

But this scheme is valid only for the perturbation whose wavelength is larger than the 
Hubble horizon scale. For the matter field which satisfies energy conditions, the perturba- 
tion terms induced by the spatial gradient terms grows in time and finally dominates the 
background solution. This occurs when the wavelength of the perturbation equals the Hub- 
ble horizon scale. After this time, the wavelength of the perturbation becomes shorter than 
the horizon and the result of the gradient expansion becomes unreliable. 

The similar situation occurs in the field of non-linear dynamical systems. To obtain 
temporal evolution of the solution of a non-linear differential equation, we usually apply 
perturbative expansion. But naive perturbation often yields secular terms due to the res- 
onance phenomena. The secular terms prevent us from getting approximate but global 
solutions. There are many techniques to circumvent the problem, for example, averaging 
method, multi-time scale method, WKB method and so on Although these method 

yield globally valid solutions, they provide no systematic procedure for general dynamical 
systems because we must select a suitable assumption on the structure of perturbation series. 

Renormalization group method [f7|-p]] as a tool for global asymptotic analysis of the 
solution to differential equations unifies the techniques listed above, and can treat many 
systems irrespective of their features. Starting from a naive perturbative expansion, the 
secular divergence is absorbed to the constants of integration contained in the zeroth order 
solution by the renormalization procedure. The renormalized constants obey the renormal- 
ization group equation. This method can be viewed as a tool of system reduction. The 
renormalization group equation corresponds to the amplitude equation which describes slow 
motion dynamics in the original system. We can describe complicated dynamics contained 
in the original equation by extracting a simpler representation using the renormalization 
group method. 

In this paper, we apply the renormalization group method to the gradient expansion 
of the Einstein equation. Our purpose is to obtain the renormalized long-wavelength solu- 
tion of the Einstein equation which is also valid for late time. Through the procedure of 
renormalization, we extract slow motion from the Einstein equation. 

This paper is organized as follows. In section 2, we briefly review the gradient expan- 
sion. In section 3, we introduce the renormalization group method by using two examples. 
In section 4, the renormalization group method is applied to the solution of the gradient 
expansion. In section 5, we solve the renormalization group equation for several situations. 
Section 6 is devoted to summary and discussion. We use the unit in which c = h = SnG = 1 
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throughout the paper. 



II. LONG- WAVELENGTH EXPANSION OF EINSTEIN+DUST SYSTEM 

We use synchronous reference frame: 

ds^ = -dt^ + -iij{t, x)dx'dx^. (1) 

The Einstein equation with dust fluid is 

+ - KfKi = 2p, (2) 
K^i - Klj = v^l + u'^Uk p Ui, (3) 



where K-ij is the extrinsic curvature, p is the energy density of dust, Ui is the three velocity 
of dust. If we neglect terms with spatial derivative, the solution can be written 

^f^a\t)h,^{x), p^'^^poit), uf^O, (6) 

where hij (x) is an arbitrary function of spatial coordinate (seed metric) , and the scale factor 
a{t) and the energy density po(t) obey usual spatially flat FRW equation. Using this as the 
zeroth order solution, we solve this system by the following gradient expansion: 

-fij = a\t)hij{x) + {f2{t)Rij{h) + g2{t)R{h)Kj) + • • • , (7) 
Ui = U3{t)ViR{h) + ■■■, (8) 
p^Po(t)+p2(t)R(h) + ---, (9) 

where Rij{h) is the Ricci tensor of the seed metric hij. The zeroth order solution is 

7l?=t^/%.(.), .f = 0, p('') = ±. (10) 



The solution to the next order is 



9tV^ 1 



7lf = -:^(^i?,,(/i)-^i?(/i)/.,,), (11) 



The solution up to the second order of spatial gradient can be written as follows: 

9 o/, /„ ... 1 



t'" {hii - f [RiiW - } , (12) 



We can see that the perturbation term, which originated from the spatial gradient of the seed 
metric, grows as the universe expands and finally has the same amplitude as the background 
term at t ~ when the wavelength of the perturbation equals Hubble horizon scale. 

After this time, the wavelength of perturbation becomes smaller than horizon scale and the 
long-wavelength expansion breaks down. To make the gradient expansion applicable to the 
perturbation whose wavelength is smaller than the horizon scale, we use the renormalization 
group method. 



III. RENORMALIZATION GROUP METHOD 

The renormalization group method improve the long time behavior of naive per- 

turbative expansion. We explain the basic concept of the renormalization group method 
using two examples. First one is a harmonic oscillator. The equation of motion is 

X + X = —ex, (13) 

where e is a small parameter. We solve this equation perturbatively by expanding the 
solution with respect to e: 

X = Xq + exi + ■ ■ ■ . (14) 

The solution up to 0(e) becomes 

X = Bocost + Cosint + - {t - to){CoCost - Bosint) + O(e^), (15) 

where Bq and Cq are constants of integration determined by the initial condition at arbitrary 
time t = tQ. This naive perturbation breaks down when e {t — to) > 1 because of the secular 
term. To regularize the perturbation series, we introduce an arbitrary time n, split t — to 
as t — fi + fi — to, and absorb the divergent term containing fi — to into the renormalized 
counterparts B and C of Bq and Cq, respectively. 
We introduce renormalized constants as follows: 

Bo = B{fi) + e6B{fi,to), Co = C{fi) + e6C{fi,to), (16) 

where 6B and 6C are counter terms that absorb the terms containing fi — to in naive solution. 
Substitute Eq. ([T6|) to Eq. (p^Sf) , we have 

X = B{fi) cost + C{fi) sint 

+e S^6B cost + 6C sint + ^{t - ^ + fi - to){C{fi) cost - B{n) sint)| . (17) 

SB and SC are chosen as 

5B{fi, to) + ^(/i - to)Cifi) = 0, 5Cifi, to) - ^(/i - to)B{fi) = 0. (18) 

Using the relation e6B = Bq — B{n), eSC = Co — C(/i), 
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B{to) = B{fx) - - to)C(/x), C(to) = C{ii) +'-{^^- t,)B{ii). (19) 

These equations define the transformation up to 0(e): 

n,^t,: (i?(to),C(to))-(i?(/i),0(/x)). 
Explicit form of the transformation is 

TZ satisfies the composition law: 

n 6^n ^ f(/^-/^i)W 1 f(/^i-^o) V 2x 

i-f(/^-/^i) 1 jl-f(-"i-^o) 1 

^( 1 K/^-^o)' 

l^-f(/i-to) 1 
= ^M-io, (21) 



0(6^) 



and TZq = 1 and 7^_^ is the inverse transformation of 7^^. So we can conclude that the 
transformation forms a Lie group up to 0(e). Assuming the properties of Lie group, we 
can extend the locally valid expression (^) to a global one, which is valid for arbitrary 
large /i — to- We apply this transformation to get {B (n) , C {fj,)) at arbitrary large /x. By 
differentiating Eq. (|19|) with respect to fj, and setting to = fj,, we have renormalization group 
equations: 

The renormalized solution becomes 

x{t) = B{^i) cost + 0(/i) sint + ^ (^ - Ai)(0(At) cost - E(/i) sint). (23) 



Solving the renormalization group equation (|2^) and equating /i and t in (|2^) eliminates the 
secular term and we get uniformly valid result: 

x{t) = 5(0) cos (l + ^ + C(0) sin (^1 + t. (24) 

Second example is the Einstein equation for the FRW universe with dust. The spatial 
component of the Einstein equation is 

a + ^(d)^ = -|e-2", (25) 

where a{t) is the logarithm of the scale factor of the universe a{t) and e is the sign of the 
spatial curvature. The exact solution is given by 
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flo (1 — COST]), t = ao{r] — sin rj) for e = 1 
a{t) = e"(*) = { ao^, t = ao^ for e = (26) 

ao (coshr] — 1), t = ao (sinh?7 — rj) for e = — 1 



We solve Eq. (|25D perturbatively by assuming that the right hand side is small. This is 
expansion with respect to small spatial curvature around the fiat universe. By substituting 
a = ao + e ai + ■ ■ ■ to Eq. (^5]), we have naive solution: 

a = lnT + Co-e ^e-2^«(r - tq) + 0{e'), (27) 

where we introduced a new time variable r = t^/^ and Co is a constant of integration 
determined by the initial condition at r = tq. 0(e) term is secular and we regularize this 
term by introducing arbitrary time fi and renormalized constant Cq = C{^) + e6C{fi, Tq): 

« = In r + C(/i) + e 5C(/i, Tq) - e ^e'^^^^^ (r - + ^ - tq) . (28) 

The counter term 6C is determined by so as to absorb ^ — Tq dependent term: 

5C(^,ro)-^e-2^(^)(^-ro) = 0. (29) 

This defines the renormalization group transformation TZ^-tq '■ C{tq) —>■ C{fi) 

C(^) = C(ro)-6Ae-2C(M)(^_^^), (30) 

and this transformation forms Lie group up to 0(e). So we can have C{fi) for arbitrary large 
value of /i — To by assuming the property of Lie group, and this makes possible to produce 
globally uniform approximated solution of the original equation. The renormalization group 
equation is 



(9/i 20 
and the solution is 



e^'^^^\ (31) 



C(/^) = ^ In (c - ^/i) , (32) 
where c is a constant of integration. The renormalized scale factor is given by 

a{t) = = T e^(^) = t^/s L - j . (33) 

As the zeroth order solution, it is possible to include another integration constant to 
which defines the origin of cosmic time t. By requiring that the renormalization group 
transformation forms Lie group, it can be shown that to does not get renormalization. So it 
is sufficient to consider the solution with the boundary condition a{t = 0) = which fixes 
the value of to equals zero. This point is different from the example of harmonic oscillator, 
in which case two integration constants B and C both get renormalization. 

We compare the renormalized solution (^) with the exact solution ( PBD and the naive 
solution (^) for the case of a closed universe (e = 1). The scale factor of the exact solution 
has maximum at t = vr and goes to zero at t = 2n. The naive solution does not show 
this behavior. The renormalized solution improves the naive solution and reproduces the 
expanding and contracting feature of the exact solution (Fig.l). 
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IV. RENORMALIZATION OF LONG- WAVELENGTH SOLUTION 



Now we proceed to the long-wavelength solution of Einstein equation. We renormalize 
the secular behavior of three metric jij. By introducing a new time variable r = t^^^, the 
longwavelength solution can be written as 

= T - ^ (i?., (/.) - \Rih)h,,'^ } . (34) 

We renormalize the secular term "^{Rij — \Rhij). We introduce the initial time Tq by re- 
defining the seed metric /i^, and define renormalized metric and the counter term hij{x) = 
hij{x,fi) + 6hij{x,fj,,To): 

9, . /„ 1 



9 / 1 \ 

hij{x) - -(r - To) i^R,j - -R{h)hijj 



9 / 1 \ 

= hij{x, jj,) + 6hij - -{t - fi + jj, - To) iRij - -R{h)hij] . (35) 

The counter term is determined so that it absorbs terms containing — tq: 

Shi.ifx, t) - lifi - To) f%(/i) - ]R{h)hij] = 0. (36) 



Using the definition of 6hij, 

9 / 1 \ 

hij{x,fi) = hij{x,To) - -{fi - To) iRij{h{fi)) - -R{h{fi))h,j{fi)] . (37) 

This equation defines the renormalization group transformation '■ hijijo) —>■ hij{fi) 

and this transformation is Lie group up to 0(e). We therefore can get the value of hij{fi) for 
arbitrary fi using the relation (|37|) by assuming the property of Lie group. The renormal- 
ization group equation is obtained by differentiating Eq. ( pTj ) with respect to fi and setting 
To = /i: 



d 9 

— /i,,(x,/i) = -- 



Rij{Kii)) - ^R{h{fi))hij{fi) 



(38) 



and renormalized solution is 



-fij = t^/^hij{x,t), (39) 
3t2 + ^ 



-^ + ^^Rihix,t)). (40) 



V. SOLUTION OF RENORMALIZATION GROUP EQUATION 

We solve the renormalization group equation (|38D for some special cases and see how the 
renormalization group method improves the behavior of the long-wavelength solution. 
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A. FRW case 



The metric is 



hij (x, t) = Q"^ (t) aij (x) , Rij (a) = ^aij R{a) , 



where (Jij{x) is the metric of three dimensional maximally symmetric space. In this case, 
the renormalization group equation reduces to 

^^^'W = -^fc, (A; = ±1,0), 

and the solution is 



where c is a constant of integration. The renormalized solution is 

7, =t^/^ (c-^-^f/']a.,ix 



± (i EJifH 

3t^[ + 20 c - f t2/3 
and the scale factor of the universe is given by 



P=^[^ + ^ . 9fc../J ^ (41) 



a(t)=t^/yc-^t2/3. (42) 

This is the same as the solution 

B. spherically symmetric case 

In spherical coordinate (r, 6',0), the metric is 

h,^ = diag {a\t, r), B\t, r), B\t, r) sin^ 0) . (43) 
The renormalization group equation (|38| ) becomes 

2i?|^ = ^(-^ + ^)- (^^ component) (45) 
The solution of this equation is given by 

B=(l- ^^t] ' /?(r), (46) 



10 ; 

B J. 

^Jl - a{r) (3^{ry 



A = , (47) 
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where a{r) and (3{r) are arbitrary function of r. The renormahzed metric and density are 



, 2 

, 2 



ds-" = -de + \ dr' + (e'^BYdnl (48) 

1 — a ^ ^ 



4 



3 


3a (l 




10 ' 


+ a,^/3 




5^4/3 


(1 - 






(1- 


10 j 




~ 20*^''' 



(49) 



By choosing = r, a = 0, ±1, we can recover the solution of FRW case (^Tl). This solution 
corresponds to Toleman-Bondi solution ||12||: 



where 



ds" = -de + ^'"^I'l), , dr^ + r)'cin^ (50) 
1 + 2hi[r) 



-^(1 - cosr^), t = (.2^^)3/2 - sinr^) for E(r) < 
i?(t,r) = <j (9Ml))'/%2/3^ fo, ^(^) ^ (51) 

^(coshr/-l), ^= (2i?ff)l/2 (sinh r^-r^) for E(r) > 



and 



, , 2MJr) 

E{r), M{r) are arbitrary function of radial coordinate r and are connected to the initial 
distribution of the spatial curvature and the initial distribution of mass density of dust, 
respectively. They have the following correspondence with the functions a,/3 contained in 
the renormalization group solution 



2E{r) ^ -a{r)(3\r), ^M(r) ^ /3^{r). (53) 

The renormalized solution well reproduces the feature of the metric of spherically symmetric 
gravitational collapse of dust. 



C. Szekeres solution 

In Cartesian coordinate {x,y, z), the metric is assumed to be 

/lij = diag(l,l,v42(r,x,?/, z)) . (54) 
The renormalization group equation (|3^) reduces to the following three equations: 

= J^, {xy component) 

A —A 

= — '^^ — (xx component) (55) 



dA^ 1 



dr 2 



--A{A^^^ + A^yy), {zz component) 



The solution of this equation is 



A = g{z){x^ + y^) + ^g{z) T + c{z), 
5 

where c, g are arbitrary function of z. The renormahzed metric is 



dx'^ + dy^+ g{z) (x^ + y^) + 



e''^ + ciz) 



dz' 



(56) 



(57) 



This is the Szekeres' exact solution |T^, which represents one dimensional gravitational col- 
lapse. It is known that the "naive" gradient expansion reproduces this solution by including 
the fourth order spatial gradient . We obtained the solution using the second order spatial 
gradient with renormalization. In this case, renormalization procedure much improves the 
naive solution. 



VI. SUMMARY AND DISCUSSION 

We have investigated renormalization of the long-wavelength solution of the Einstein 
equation which has secular behavior. We applied the renormalization group method to 
improve the long time behavior of the solution. We introduced renormalization point fi ~ 
L/H^^ where L is physical wavelength of perturbation. For /i ^ 1(L ^ H^^), the naive 
expansion (gradient expansion) well approximate the exact solution. But as the universe 
expands, jj, ~ / L oc t^^^ becomes large and at /x ~ 1(L ~ H~^), perturbative expansion 
breaks down. After this time, the wavelength of the perturbation is smaller than the horizon 
scale. Using renormalization transformation which has the property of Lie group, we can 
renormalize the secular terms caused by the spatial gradient of the background seed metric. 
The renormalized metric can be extended to large value of /i, which corresponds to the small 
scale fluctuation. 

We obtained the solution of the renormalization group equation for FRW, spherically 
symmetric and Szekeres cases. The behavior of the renormalized solution indicates that they 
describe the collapsing phase of the system qualitatively well. The renormalization group 
method is regarded as the procedure of system reduction. This means the renormalization 
group Eq. (|38|) is reduced version of the original Einstein equation and describes slow motion 
dynamics of the original equation. We expect the renormalization group equation ( PBD has 
physically interesting properties and solutions which contained in the Einstein equation. 

We can look the renormalization of the long-wavelength solution from the view point of 
the backreaction problem in cosmology. The naive solution represents the evolution of the 
perturbation with the fixed background metric. By renormalizing the naive solution, the 
constants htj (x) contained in the background solution becomes to have the time dependence 
due to the spatial inhomogeneity. So we can investigate how the spatial inhomogeneity affects 
"background" metric hij by solving the renormalization group equation. The remarkable 
feature of the renormalization group approach is that it does not need any spatial averaging 
which is necessary in conventional approach to the backreaction problem in cosmology [P 



In this paper, we considered only dust fluid as the matter field. But the gradient expan- 
sion can treat the perfect fluid with pressure and scalar fields. For the scalar field system, 
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although we cannot write down exact solutions of the naive expansion in general, it is possi- 
ble to define the renormalization group transformation and obtain the renormalized solution 
formally. It is interesting to apply the renormalization group method to the infiationary 
model and investigate the effect of the stochastic quantum noise on the background geom- 
etry. We can construct the inhomogeneous model of stochastic inflation pr5|-[T^ and this is 
our next subject to be explored. 
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FIGURES 
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FIG. 1. The evolution of the scale factor for a closed FRW universe with dust. The solid 

curve is the exact solution, thin dashed curve is the naive solution and thick dashed curve is the 

renormalized solution. 



13 



